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Preface 


The  primary  purpose  of  this  paper  was  to  modify  a::  air- 
to-air  missile  model  used  by  the  Cruise  Missile  Independent 
Survivability  Team  (CMIST)  of  the  Aeronautical  Systems 
Division  (ASD) .  The  model  did  not  have  incorporated  within 
it,  data  pertinent  to  the  cruise  missile.  The  modification, 
integrated  the  required  data  into  the  model  so  that  the 
model  could  be  used  to  evaluate  an  air-to-air  missile's 
performance  against  the  cruise  missile. 

The  model  itself  utilizes  a  hypercomplex  number  algebra 
which  avoids  usage  of  direction  cosine  matrices  and  enables 
real  time  of  flight  output.  Prior  to  the  completion  of  this 
thesis,  no  adequate  documentation  existed  that  explained  the 
algebra  used  and  its  translation  into  computer  code.  Bits 
and  pieces  were  gathered  together  and  combined  into  the 
appendices  of  this  work,  which  should  prove  helpful  in 
understanding  the  model. 

I  wish  to  thank  Captain  Aaron  DeWispelare  for  his 
encouragement  and  help  in  completing  this  work.  I  also  wish 
to  thank  Molly  Bustard,  my  favorite  librarian,  who  consis¬ 
tently  provided  a  bright  spot  during  many  a  dreary  day  these 
last  eighteen  months.  I  must  also  acknowledge  Bob  and  Diane 
Turelli,  who  provided  an  escape  from  the  academic  grind  when 
it  was  needed,  and  finally  my  mother,  sister  and  her  family 
who  gave  me  the  moral  help  I  needed  to  complete  my  studies. 
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Abstract 

The  air-to-air  missile  model  used  by  the  Cruise  Missile 
Independent  Survivability  Team  did  not  contain  data  useful 
to  the  cruise  missile.  The  objective  of  this  study  was  to 
modify  the  subroutine  SIGNAL  in  the  model  to  incorporate 
the  cruise  missile  data.  The  modification  required  an 
understanding  of  quaternion  algebra  utilized  within  the 
model  to  represent  three-dimensional  motion.  These  quater¬ 
nions  allow  real  time  outputs  from  the  model  for  use  by 
tactical  ranges.  The  study  contains  a  discussion  of  qua¬ 
ternions  and  their  algebra. 
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I  Introduction 


Problem  Statement 

The  Cruise  Missile  Independent  Survivability  Team 
(CMIST)  obtained  the  computer  model  of  the  air-to-air  mis¬ 
sile  performance  package  utilized  in  the  ACEVAL/AIMVAL  pro¬ 
gram.  The  model  involves  a  new  technique  of  representing 
the  mechanics  of  the  air-to-air  missile's  flight  by  using 
quaternions . 

The  hypercomplex  quaternions  were  not  used  regularly 
in  a  computer  at  the  Aeronautical  Systems  Division.  Major 
Ken  Madsen  (CMIST)  expressed  a  desire  for  an  analysis  of 
how  quaternions  are  utilized  in  the  model  to  be  able  to  make 
future  changes  if  necessary.  In  addition,  he  requested  a 
modification  to  the  infrared  signature  subroutine  (SIGNAL) 
which  would  incorporate  infrared  data  generated  for  the 
cruise  missile  and  not  currently  in  the  model. 

Background 

In  1843  Sir  William  Rowan  Hamilton  developed  a  new 
algebra  of  quadruples  of  numbers  which  he  named  quaternions. 
His  concept  was  intended  to  represent  motion  in  three-space, 
predating  the  familiar  vector  methods  of  J.  Willard  Gibbs. 
Hamilton's  quaternion  theory,  although  the  foundation  of 
vector  algebra,  was  from  the  viewpoint  of  the  student  of 
mathematical  physics,  a  confusing  mixture  of  scalars  and 
vectors.  The  vector  which  represented  three-dimensional 
motion  seemed  to  play  a  servile  role  as  part  of  the  quatern¬ 
ion.  The  real  power  of  the  quaternion,  however,  came  from 
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the  vector  aspects  of  the  algebra  which  was  later  extracted 
by  J.  Willard  Gibbs.  Gibbs  simplified  the  quaternion  meth¬ 
ods  of  Hamilton  and  originated  the  vector  analysis  techniques 
that  we  are  familiar  with  today  (Ref  3) . 

Because  the  concept. of  a  three-dimensional  vector  was 
more  acceptable  than  the  hyper complex  quaternion  (Ref  9) , 
vector  algebra,  a  subset  of  quaternion  algebra,  was  adopted 
by  the  mathematical  community  to  describe  three-dimensional 
motion.  With  the  advent  of  recent  computer  missile  simula¬ 
tions,  however,  different  computational  methods  were 
required  to  reduce  the  computation  times  to  the  real  times 
of  flight.  To  accomplish  this,  the  concept  of  quaternion 
algebra  was  resurrected  and  found  capable  of  not  only  reduc¬ 
ing  computation  time  but  also  storage  space  normally  required 
by  the  nine  element  direction  cosine  matrices. 
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II  Subroutine  SIGNAL 


The  missile  model  obtained  by  CMIST  was  designed  to 
work  with  actual  launch  and  target  aircraft  in  mock  combat 
on  a  controlled  range.  While  the  aircraft  maneuver,  data 
is  relayed  to  a  range  computer  by  telemetry  pods.  The  com¬ 
puter  tracks  and  stores  their  flight  profiles  and  also  does 
real  time  missile  launch  calculations  using  the  missile 
model  to  determine  if  a  hit  or  miss  is  achieved. 

The  missile  model  acquired  by  CMIST  was  to  do  some 
effectiveness  studies  of  several  air-to-air  missiles  against 
the  cruise  missile.  The  model,  however,  did  not  contain 
the  required  data  to  calculate  the  infrared  signal  generated 
by  the  cruise  missile  and  could  not  be  used  until  the  data 
was  incorporated  into  the  model. 

The  infrared  signal  is  calculated  within  the  model  in 
a  subroutine  called  SIGNAL.  The  subroutine  then  uses  the 
approximation  assumption  of  the  basic  model  that  the  signal 
strength  is  directly  proportional  to  the  maximum  line  of 
sight  tracking  rate.  Utilizing  this  assumption,  the  sub¬ 
routine  compares  the  calculated  signal  to  several  threshold 
values  to  determine  if  the  strength  of  the  signal  is  adequate 
for  the  gimbaled  seeker  head  to  track  the  target  aircraft  at 
the  calculated  line  of  sight  rate. 

The  modification  proposed  not  only  involved  the  incor¬ 
poration  of  the  infrared  data,  but  also  required  a  redefi¬ 
nition  of  the  line  of  sight  reference.  This  involved  deter¬ 
mining  the  line  of  sight  in  the  target  coordinate  system  by 

means  discussed  in  the  Methodology. 
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Ill  Methodology 


Two  major  problems  were  encountered  in  the  modification 
of  the  CMIST  missile  model.  First,  the  model  did  not  uti¬ 
lize  standard  direction  cosine  coordinate  transformations, 
an  accepted  technique  to  represent  spatial  relationships 
among  several  coordinate  systems,  in  SIGNAL.  Secondly,  the 
method  employed  by  SIGNAL  involved  hypercomplex  quaternions 
which  are  foreign  to  many  engineering  students.  (The 
interested  reader  will  find  a  discussion  of  quaternion 
development  and  their  properties  in  Appendix  A.)  This  sec¬ 
tion  discusses  the  relative  merits  of  direction  cosines  and 
quaternions,  with  a  justification  for  the  methodology  selec¬ 
ted. 

Direction  Cosines 

Direction  cosines  have  been  the  accepted  engineering 
technique  to  represent  spatial  relationships  among  several 
coordinate  systems.  They  are  relatively  easy  to  use 
because  direction  cosines  follow  the  rules  of  conventional 
vector  algebra  and  can  be  stored  in  a  computer  as  elements 
of  an  array.  A  vector  can  then  be  defined  in  another  coor¬ 
dinate  system  by  utilizing  the  array  as  a  directional  cosine 
matrix  and  using  matrix  multiplication  to  transform  the  vec¬ 
tor  coordinates.  For  example,  the  vector  x  will  be  redefined 
in  some  primed  coordinate  system  by  multiplying  it  by  the 
Direction  Cosine  Matrix  (DCM)  which  links  the  unprimed  coor¬ 
dinate  system  to  the  primed  system 
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X<DCM) 


*  x' 


Direction  cosines,  although  familiar  and  relatively 
easy  to  use  on  a  computer  also  have  several  disadvantages. 
These  consist  of 

1)  Each  matrix  requires  a  nine  element  storage  allo¬ 
cation. 

2)  Trigonometric  functions  are  required  to  form  the 
elements . 

3)  Many  multiplications  of  the  elements  must  be  per¬ 
formed  in  a  proper  order. 

The  direction  cosine  method  would  result  in  a  major  rewrite 
of  the  subroutine  SIGNAL  but  would  enable  a  familiar  engi¬ 
neering  method  to  be  employed. 

Quaternion  Method 

The  quaternion  method  employed  by  the  model  derives  a 
quaternion  to  relate  one  coordinate  system  to  another.  This 
methodology  is  discussed  in  Appendix  B.  The  quaternion  q 
then  operates  on  a  vector  x  converting  it  to  x1  by 

q”1  X  q 

where  q  ^  is  the  inverse  of  q.  (The  interested  reader  is 
again  referred  to  Appendix  A.) 

Quaternions  consist  of  only  four  elements  and  thus  are 
advantageous  over  direction  cosines  since  they  require  fewer 
arithmetic  operations  and  less  storage  space.  The  disadvan¬ 
tage  of  quaternions  is  the  requirement  of  understanding  their 

unique  properties  and  operating  rules. 
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Selection 


To  avoid  a  major  revision  of  the  computer  code,  quater¬ 
nions  were  selected  to  modify  SIGNAL.  First,  quaternions 
did  not  require  a  conversion  of  the  vectors  from  the  qua- 
ternion  factored  format  into  the  conventional  format  used  by 
direction  cosines.  Secondly,  after  the  operations  were  com¬ 
pleted  the  resulting  vectors  would  have  to  be  reconverted  to 
the  factored  format  used  by  the  rest  of  the  program.  Fin¬ 
ally,  the  construction  of  direction  cosine  matrices  could 
be  avoided  along  with  a  dimensioning  of  their  storage  space. 
Quaternions  were  perceived  to  be  the  way  to  go  because  their 
methods  would  have  to  be  understood  in  either  case. 
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T.V  Modification  of  SIGNAL 


Data  Incorporation 

The  modification  of  the  subroutine  was  to  incorporate 
the  infrared  data  supplied  by  CMIST.  The  data  used  was 
generated  by  a  computer  program  into  a  set  of  data  points 
at  various  azimuths  and  elevations  about  the  cruise  missile. 
This  data  was  then  incorporated  into  several  arrays  which 
served  as  a  tabulated  reference ,  the  idea  being  to  enter  an 
array  with  an  appropriate  range  and  angle  off  from  the  tar¬ 
get  and  be  able  to  select  the  relative  infrared  signal 
strength. 

Coding 

The  principal  problem  in  obtaining  the  infrared  signal 
was  to  convert  the  missile  seeker  line  of  sight  to  the  tar¬ 
get  coordinate  system  so  the  tabulated  data  could  be  used. 
The  range  vector  is  defined  in  the  universal  coordinate 
system.  This  same  vector  in  the  quaternion  factored  format 
is: 


r  *  (rx  +  kZ)  i  . 

A  quaternion  q,  which  relates  the  universal  coordinate  sys¬ 
tem  to  that  of  the  target  coordinate  system  is  used  to 
redefine  r  in  the  target  system  by 

q”1  r  q  . 

The  interested  reader  is  referred  to  the  methods  of  Appendix  3. 
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With  the  range  vector  defined  in  target  coordinates,  the 
relative  position  of  the  seeker  line  of  sight  in  azimuth  and 
elevation  can  be  determined.  With  this  relative  position, 
a  data  point  can  be  computed  from  the  arrayed  information. 
Using  this  data  point,  the  received  signal  strength  is 
reduced  by  a  factor  of  one  over  the  range  squared.  Upon 
determination  of  the  strength  of  the  received  signal,  the 
program  reverts  to  the  original  line  of  sight  rate  tolerance 
calculations  and  the  simulation  continues. 

Validation 

The  modification  to  SIGNAL  was  validated  by  comparing 
the  infrared  signature  output  of  the  old  model  with  that  of 
the  modified  routine.  To  accomplish  this,  the  function 
which  formerly  computed  the  infrared  signal  based  upon 
aspect  angle  and  range  was  extracted  from  the  original 
model.  Using  this  extracted  function,  a  series  of  data 
points  were  generated  to  fill  the  azimuth  and  elevation  data 
arrays  used  by  the  modified  routine.  Once  the  arrays  were 
filled,  infrared  signal  strengths  were  generated  using  both 
models  for  various  aspects  and  ranges  of  three,  six  and 
nine  thousand  feet. 

Old  Model 

The  original  infrared  data  generating  function  used  in 
the  old  model,  utilized  data  generated  for  a  particular  air¬ 
craft  by  the  Navy  at  China  Lake.  This  data  was  plotted  into 
several  cardioid  shaped  curves  (Figure  1)  which  varied  in 
size  as  a  function  of  range  and  altitude.  A  function  was 
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Fig.  1.  IR  Signature  Cardioid  (Watts/Steradian) 


fit  to  these  curves  for  different  aspects  with  parameters 
of  atmospheric  attenuation,  range  and  the  cosine  of  the 
aspect  angle.  The  old  model  used  this  functional  format  to 
determine  the  infrared  signal.  A  unique  feature  of  this 
function  was  that  the  output  was  in  the  form  of  decibels 
which  are  called  Phasey  dB's.  The  data  supplied  by  China 
Lake  was  in  decibels  formed  by  taking  the  natural  logarithm 
of  the  signal  and  multiplying  it  by  a  factor  of  ten  rather 
than  the  standard  method  of  ten  times  the  logarithm  (base  10) 
of  the  signal.  The  model  uses  these  Phasey  dB's  in  all  of 
its  signal  comparisons,  something  to  be  aware  of  when  try¬ 
ing  to  follow  the  logic. 

Comparing  Results 

Once  the  old  infrared  function  was  understood,  data 
points  were  generated  for  1000,  5000  and  13,000  feet  ranges 
at  various  aspects  and  elevations.  These  data  points  (.in 
watts/steradian)  were  next  inserted  into  the  azimuth  and 
elevation  arrays  utilized  by  the  modified  routine.  (The 
interested  reader  is  referred  to  Appendix  D . ) 

Once  both  infrared  routines  were  set  up  using  informa¬ 
tion  based  upon  the  same  aircraft,  a  comparison  set  of  data 
was  generated  by  both  methods  to  determine  the  accuracy  of 
the  new  infrared  routine.  The  results  are  as  shown  in 
Table  I. 

Some  typical  launch  ranges  for  am  IR  missile  at  various 
aspects  in  the  rear  hemisphere  of  the  target  (.180°  being  a 
direct  tail  shot)  were  selected.  A  comparison  of  the  old 
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and  new  IR  signature  routines  shows  that  the  new  model  var¬ 
ies  at  most  by  less  than  two-and-one-half  percent  from  the 
old  and  that  at  an  extreme  aspect  angle. 


TABLE  I 

Comparison  of  Old  and  New  Models 


IR  Signal 

(Phasey  dB)  % 


Range 

Aspect 

Old 

New 

Difference 

3,000 

180 

102.8 

103.4 

0.6 

6,000 

180 

87.9 

88.0 

0.1 

9,000 

180 

79.1 

79.4 

0.3 

3,000 

170 

104.2 

103.6 

0.6 

6,000 

170 

89.0 

87.7 

1.5 

9,000 

170 

80.0 

78.9 

1.4 

3,000 

165 

102.7 

103.7 

1.0 

6,000 

165 

87.3 

87.5 

0.2 

9,000 

165 

78.3 

78.7 

0.5 

3,000 

150 

98.0 

99.2 

1.2 

6,000 

150 

82.2 

82.5 

0.4 

9,000 

150 

73.0 

73.5 

0.7 

3,000 

140 

94.6 

95.8 

1.3 

6,000 

140 

78.7 

78.7 

0.0 

9,000 

140 

69.3 

69.6 

0.4 

3,000 

135 

92.9 

94.3 

1.5 

6,000 

135 

76.8 

77.1 

0.4 

9,000 

135 

67.3 

67.9 

0.9 

3,000 

120 

87.2 

88.8 

1.8 

6,000 

120 

70.7 

71.1 

0.6 

9,000 

120 

61.0 

61.6 

1.0 

105 

83.1 

84.9 

2.2 

105 

66.3 

66.7 

0.6 

9,000 

105 

56.4 

57.1 

1.2 

3,000 

90 

82.6 

84.6 

2.4 

90 

65.5 

65.9 

0.6 

9,000 

90 

55.4 

56.1 

1.3 
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V  Summary  and  Conclusions 


Quaternions,  although  initially  appearing  to  be  a  com¬ 
bination  of  a  vector  and  a  scalar,  possess  several  desirable 
properties  when  compared  to  conventional  vector  methods. 

These  properties  have  been  applied  in  a  missile  model  which 
has  been  acquired  by  the  CMIST  branch  at  ASD.  The  model 
was  obtained  to  evaluate  several  air-to-air  missiles  against 
a  cruise  missile  target. 

To  use  the  model,  data  which  represent  the  cruise  mis¬ 
sile  characteristics  had  to  be  incorporated.  This  required 
a  working  knowledge  of  some  of  the  internal  mechanics  of  the 
model  and  how  the  quaternions  were  utilized.  A  working  know¬ 
ledge  of  the  model's  architecture  and  coding  was  also  required 
to  modify  the  subroutine  SIGNAL  which  would  take  infrared  data 
supplied  by  CMIST  and  integrate  it  with  the  existing  model. 
This  research  accomplished  the  incorporation  of  the  required 
data  and  compared  the  results  with  the  infrared  data  gener¬ 
ating  function  which  the  model  has  originally  incorporated. 

The  most  significant  result  of  this  research  was  the 
modification  to  the  subroutine  SIGNAL  and  a  greater 
detailed  description  of  how  quaternions  are  used  in  the  model 
than  available  in  the  documentation.  The  modification  of 
SIGNAL  saved  the  Air  Force  from  contracting  the  work  out  and 
is  currently  capable  of  generating  the  required  data  to  com¬ 
plete  an  evaluation  of  the  cruise  missile.  Also,  utilizing 
the  detailed  description  of  the  computer  code  along  with  the 
description  of  how  it  applies  to  the  quaternions,  further 
modifications  should  be  possible  if  necessary  in  the  future. 
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Appendix  A 
Quaternions 


Initial  Development 

Quaternions  are  specified  using  four  real  parameters, 
a  scalar  and  three  other  units  i,  j  and  k.  The  four  param¬ 
eters  specify  the  properties  that  a  quaternion  possesses 
when  it  acts  as  an  operator.  These  properties  allow  for  the 
rotation  of  a  line  vector  through  a  given  angle  about  an  axis 
through  its  origin  and  a  stretching  of  the  vector  by  a  given 
ratio.  Of  the  four  parameters;  two  are  required  to  specify 
the  axis  of  rotation,  one  to  specify  the  angle  of  rotation 
and  one  to  specify  the  ratio  of  stretch  (Ref  8) . 

A  quaternion  q  can  be  expressed  in  the  form 

q  =  q0  +  qLi  +  q2j  +  q3*  • 

Thus  it  consists  of  two  parts,  a  scalar  qg  and  a  vector 

q  =  qxi  +  q2j  +  <33k 

where  q^,  q2  and  q3  are  the  rectangular  Cartesian  coordinates 
of  a  point  P  and  i ,  j ,  k  represent  unit  vectors  in  the  posi¬ 
tive  direction  of  the  x,  y  and  z  axes,  respectively.  The  vec¬ 
tor  defines  a  line-vector  from  the  origin  0  to  the  point  P 
(Ref  8) . 

Two  quaternions  q  and  q'  are  defined  to  be  equal  if  (Ref  2) 

q  -  q0  +  +  q2j  +  q3* 

q'  *  qj  +  qji  +  q2j  +  q3* 
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and 


q0  '  qi  =  »  ^2  *  q2  and  q3  =  q3 


q0  = 

There  also  exists  a  zero  quaternion  {Ref  2) 

0  +  Oi  +  Oj  +  Ok  or  simply  (0,0, 0,0) 

such  that 

q  +  10,0,0,0)  =■  q0  +  qji  +  q2j  *  -  q  . 

Two  quaternions  are  added  by  adding  their  scalar  parts 
and  the  corresponding  coefficients  of  their  vector  parts 
(Ref  2) 

q  +  q*  =  (q0+qQ)  +  (qi+qi)i  +  Cq2+q2)j  +  (q3+q3)k  • 

Thus  the  sum  of  two  quaternions  is  also  a  quaternion. 

For  any  scalar  X  a  quaternion  cam  be  defined  (Ref  2) 

Xq  =  X(qQ  +  q1i  +  q£j  +  q3k) 

~  xqQ  +  ^qxi  +  Xq2j  +  Xq3k  . 

Also  for  X  =  -1 

(-l)q  =  -qQ  -  qxi  -  q2j  -  q3k 

and 

q  +  (-l)q'  *  (-<3o“qO)  +  tql"qi)jL  +  Cq2_q2)j  +  Cq3_q3>k  . 

So,  quaternions  obey  the  algebraic  laws  of  addition  and 
scalar  multiplication  (Ref  2) . 
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It  can  also  be  shown  that  for  any  quaternions  p ,  q  and 
r  and  any  scalars  y  and  X  (Ref  2) 

p  +  q  =  q  +  P 

(p+q)  +  r  =  p  +  Cq+r) 

Xq  =  qX  ,  (Xy )  q  =  X(yq) 

(X+y ) q  =  Xq  +  yq 

and  X  (p+q)  =  Xp  +  Xq  . 


To  examine  quaternion  multiplication,  the  vector  ele¬ 
ments  of  a  quaternion  must  be  examined.  In  a  Cartesian 
coordinate  system,  three  mutually  perpendicular  axes  exist: 
x,  y  and  2.  For  a  right-handed  system,  the  positive  direc¬ 
tions  are  as  shown  in  Figure  2.  Rotation  is  defined  as 
positive  when  viewed  from  the  origin  (0)  as  a  clockwise 
rotation  of  the  coordinate  axis.  Thus  rotation  from  x  to  v 
about  z,  y  to  z  about  x,  and  z  to  x  about  y  are  all  positive 


Fig.  2.  Cartesian  Coordinate  System 
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The  vectors  i,  j  and  k  are  defined  as  unit  vectors 
oriented  along  the  x,  y  and  z  axes,  with  each  specifying  a 
direction  along  its  respective  axis  and  possessing  a  magni¬ 
tude  of  one.  The  rotation  of  the  axis  can  be  defined  by 
the  multiplication  of  two  of  the  unit  vectors,  i.e.  the 
multiplication  of  i  into  j  or  ij  is  defined  to  be  the  turn¬ 
ing  of  j  by  +90°  in  the  plane  normal  to  i  (Ref  6) .  There¬ 
fore 


Similarly, 


ij  =  k  . 


jk  =  i  ,  ki  =  j 


and  ji  =  -k  ,  kj  =  -i  ,  ik  =  -j  . 

Furthermore ,  since 


ij 

=  k 

and 

ik 

=  -j 

then 

-j  -  ik  - 

i  Cij  1  =  i2j 

or 

i2j 

=  -j 

therefore 

i2  - 

=  -1  . 

Similarly 

j2  =  k2  = 

ijk  =  -1  . 

All  the  familiar  algebraic  rules  of  multiplication 
except  for  the  commutative  law  are  valid  for  quaternions. 
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This  relates  to  the  fact  that  the  order  in  forming  the  pro¬ 
ducts  of  i,  j  and  k  determines  positive  and  negative  rota¬ 
tions  and,  therefore,  the  order  must  be  preserved  (Ref  6). 

Prom  the  previous  discussion  a  quaternion  product  can 
be  defined  (Ref  2) . 

qr  =  (qQ  +  q^  +  q2j  +  q3k)  (rQ  +  ^i  +  r2j  +  r-jk) 

"  q0r0  "  qlrl  “  q2r2  "  q3r3  q0  (rli  +  r2 j  +  r3k) 

+  r0(q;Li  +  q2j  4-  q3k)  +  (q2r3  -  q3r2)i 

+  ^3ri  “  ^ir3^  +  ^lr2  ~  q2rl)k 

=  q0r0  "  qlrl  "  q2r2  “  q3r3  +  ^O^l1  +  r2^  +  r3k) 

i  j  k 

+  rQ  tqii  +  q2 j  +  q3k)  +  q2  q3  . 

rl  r2  r3 

Now  consider  a  vector  v,  which  is  a  special  case  of  a  qua¬ 
ternion  whose  scalar  portion  is  zero 

v  =  0  +  ai  +  bj  +  ck 

The  quaternion  product  of  two  vectors  v  and  v1  would  be 

vv*  =  0  -  aa'  -  bb'  -  cc'  +  0(a’i  +  b’j  +  c'k) 

i  j  k 

+  0 (ai  +  bj  +  ck)  +  a  b  c 

a'  b'  c' 
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or 


W 


{aa'  +  bb'  +  cc1) (-1)  + 


1 

a 


j  k 
b  c 


a’  b*  c 


There  exists  a  unit  sphere  with  origin  0  and  an  orthagonal 
coordinate  system  as  shown  in  Figure  3.  The  vector  v 
is  defined  as  oriented  along  i  with  vector  v1  located  in  the 
i,  j  plane,  both  with  origins  at  0  and  v'  forming  an  angle 
0  with  v.  (0  <  90°) 


Fig.  3.  Unit  Sphere 
The  vectors  are  now  defined 

v  =  ai 

and  v*  =  a'i  +  b* j  .  a '  =  a 

By  definition 

|v|  =  a  ,  |v‘  |  =  [Ca')2  +  Cb')231/2 


Sin@ 


b* _  =  b* 

[Ca'i2  +  (b')2]1/2  | v 1 
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_ _ 


cose 


[(a')2  +  (b')2]1/2  |  v 1  | 


v  s  I v 1 i  and  v*  =  |v'|cos8i  +  v’ 


A  vector  product  can  now  be  defined  as 


w'  =  |  v |  i  [|v*|(icos6  +  jsin6)] 


v  I v * | ( —  cos6  +  ksin8) 


or 


w'  =  -  I  v|  I  v'  loose  +  IvMv'lksinO 


Substituting 


w'  *  -  a(V(a,)2+(b,)2)cos6  +  a^/(a' )  2+(b* ) 2 


but 


V(a' ) 2+(b' ) 2  cose 


and 


V(a’)2+(b’)2  sine  *  b' 


so 


w*  =  -  aa'  +  ab'k  . 


1  sine j 


^ksin6 


Now  consider  the  quaternion  vector  product 


j  k 

b  c 
b'  c' 

but 

b  =  c  =  c '  =  0 


i 

vv*  =  - (aa1  +  bb '  +  cc' )  +  a 

a' 


so 

i  j  k 
a  0  0 

a*  b'  0 

or 


vv*  =  -  aa'  +  ab'k  (2) 

which  is  the  same  as  (1) .  The  product  of  two  vectors  can 
be  defined  as  a  quaternion  whose  scalar  part  is 

S (vv 1 )  =  -  aa' 

or  equivalently  -j v| | v' |cos0/  and  vector  part  is 

V (vv' )  =  ab'k 

or  equivalently  j v| j v' j ksin9 .  These  correspond  with  the 
scalar  and  vector  products  as  defined  in  the  vector  algebra 
of  J.  Willard  Gibbs  (Ref  2) .  Gibbs  defined  the  scalar  or 
dot  product  of  two  vectors  u  and  v  as 

u  .  v  =*  }u|  |  v{ COS0 

and  the  vector  or  cross  product  as 

u  X  v  =  |u| |v|sin0e 
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where  e  is  a  unit  vector  perpendicular  to  u  and  v  and  0  is 
the  angle  between  u  and  v.  Thus  a  quaternion  product  may 
be  defined  as 


qq*  =  q  X  q’  -  q  •  q'  . 

Properties  of  Quaternions 

The  conjugate  of  a  quaternion  is  similar  to  the  conju¬ 
gate  of  a  complex  number  (Ref  2) .  If 

q  =  qQ  +  qj^i  +  q2j  +  q3k 

then  the  conjugate  of  q  or  q* 


q*  =  q0  -  qxi  -  q2j  -  q3*  • 

The  product  of  a  quaternion  q  and  its  conjugate  q*  is 
defined  as  the  Norm  of  q  or  N(q) 

N(q)  =  qq*  =  (qQ+q1i+q2 j+q3k)  (qQ^i^ j-q3k) 

2  2,2  2 

"  q0  +  qx  +  q2  +  q3 

and  further 

q*q  =  qq*  (Ref  2) 


The  Norm  of  a  quaternion  is  a  scalar  and  N(q)  =0  implies 
that  qQ  =  q^  =  q2  =  q3  =  0.  If  N(q)  *  1,  q  is  called  a  unit 


quaternion  (Ref  2) . 


Prom  the  quaternion  product  of  two  vectors 

i  j  k 

vv '  =  »(aa'  +  bb’  +  ccf)  +  a  b  c 

a'  b'  c 
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the  conjugate  of  the  quaternion  formed  by  the  product  is 

i  j  k 

(vv')*  »  -(aa'  +  bb'  +  cc')  -  a  b  c 

a'  b'  c* 

or  equivalently 

i  j  k 

(vv' )*  =  v'v  =  -(aa'  +  bb'  +  cc')  +  a'  b'  c' 

a  b  c 

For  any  two  quaternions 

q  ■  (qQ  +  q)  /  q'  =  (qj  +  q') 
whose  product  is 

qq'  =  (q0  +  q)  (qj  +  q'>  =  q0qQ  +  qo$'  +  qoq  +  qq' 

the  conjugate  of  their  product  will  be 

(qq')*  =  q„q£  -  q„q'  -  q£q  +  (qq' ) * 

"  lo'Jo  •  “Jo3'  •  qoq  +  q,q  ' 

or 

(qq1)*  -  (q'  -  q’) (q0  -  q)  . 

Therefore,  the  conjugate  of  the  product  of  two  quaternions 
is  equal  to  the  product  of  the  individual  conjugates  taken 
in  reverse  order  (Ref  2) . 

(qq')*  =  (q')*(q)* 
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Prom  this  property,  the  Norm  of  the  product  of  two  qua¬ 


ternions  can  be  determined. 


N (qq ' )  =  (qq1 ) (qq' ) * 


*  q(q’ ) (q1 ) *q* 

*  qN (q ' ) q* 


Since  N (q)  is  a  scalar 


qN (q  * ) q*  =  qq*N(q’}  =  N(q)N(q') 


If  q  is  not  the  zero  quaternion  then  N (q)  ^  0 .  Using 
this  fact  the  inverse  of  q  can  now  be  defined  (Ref  6) . 


or 


and 


q* 

N(q) 


NJ2L  =  x 

N(q) 


It  can  also  be  shown  that 


Nfq"1) 


1 

N  (q)  • 


Using  the  concept  of  the  inverse  of  a  quaternion,  the  prop¬ 
erty  of  division  can  now  be  addressed.  For  any  three  qua¬ 
ternions  p,  q  and  r  where 
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qr  »  p  or  rq  =  p 

solutions  for  r  may  be  obtained  by  multiplying  both  sides 
of  the  equations  by  the  inverse  of  q.  In  this  way  two  solu¬ 
tions  are  apparent 

q-1qr  =  q_1p  ;  rqq-1  =  pq"1 

or 

=  q_1p  ;  r2  =  pq-1  . 

where  r^  is  called  the  right  hand  quotient  of  p  by  q  and 
r2  is  the  left  hand  quotient  of  p  by  q.  These  solutions  in 
general  are  different  and  demonstrate  again  that  order  must 
be  preserved  in  quaternion  multiplication  (Ref  2) . 

Finally,  from  the  Norm  of  q 

N(q)  =  qq* 

=  qo  +  qj  +  +  ^3 

=  4  *  lq|2 

and  the  magnitude  of  q  is  defined 

|q|  =  [N (q) ] 1/2  »  [qg  +  q2  +  q2  +  q2j 1/2 

I  q I  -  [qj  +  I  q|  21 1/2  •  (Ref  2) 
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Product  of  a  Quaternion  and  a  Vector 

To  review/  the  quaternion  product  of  two  vectors  u  and 
v 

* 

uv  =  U  X  V  -  U  *  V 

and  the  product  of  two  quaternions  q  and  r 

qr  =  (qQ+q)  (r0+?)  *  °3orO  +  q0*  +  r0q  +  qX*  “  q*?  * 

The  conjugate  of  q  is 

q*  =  qQ  -  q 

and  the  Norm  of  q  is 

N(q)  =  qq*  =  q2  +  |q|2 

where 

q  =  qii  +  q2j  +  q3k 

and 

|q|  =  (qj  +  +  qf) 1/2  . 

The  magnitude  of  the  quaternion  q  is 

|q|  =  (n(q>)1/2  =  (qg  +  |q|2)1/2  . 

For  a  vector  x  and  a  quaternion  q,  where 

q*  q0  q 

q  =  q0  +  q  '  q  "  NT5T  "  TTfq T 
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the  quaternion  product  of  q  and  x  is 


qx  =  qQx  +  (qXx)  -  (q*x) 
xq  =  qQx  +  (xXq)  -  <x*q) 


-1 


The  right  hand  quotient  of  q  and  x  is 

& 


q”1* 


cr*  -* 

Nt^X 


(q0  -  q)x 


NTqT 


or 


-1- 

q  x  = 


qQx  -  (qXx)  +  (q*x) 

Kjq) 


qQx  +  (xXq)  +  (q*x) 


N(q) 

The  left  hand  quotient  of  x  and  q”1  is 

x(qQ  -  q) 


-1 

xq 


w 


qQx  -  (xXq)  +  (x*q) 


or 


-i  qox  +  +  ^x*q^ 

xq~  =  - 


NTqT 


In  general,  the  product  of  a  vector  and  a  quaternion  yields 
a  quaternion  whose  scalar  part  is  the  dot  product  (x*q)  and 
whose  vector  part  is  the  scalar  of  the  quaternion  times  the 
vector  and  the  vector  cross  product  (qQx  +  (qXx) ) .  However, 
in  the  case  where  q  and  x  are  perpendicular  (q*x  =  0) ,  and 
the  vector  qgX  +  qXx  is  obtained.  In  this  case  both  q,„x 
and  qXx  sure  normal  to  each  other  and  the  resultant  vector 
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qQx  +  qXx  is  normal  to  q.  Figure  4  shows  the  geometrical 
relationships  of  xq  and  q  oriented  into  the  page  (Ref  5)  . 


It  can  be  seen  that 

tana 


also 


Ixllgl  _  13L 
q0|S|  9o 


xq|  - 


t(|x|  |q|)2  +  Cq0|S|)2]1/2 
t  |x|  2|q|  2  +  q2|S|2]1/2 
[|5|2Cq2  +  |q|2)]1/2 
t|S|2|q|2]1/2 


xq| 


*1 1  q  I 


Figure  4  shows  that  a  is  a  negative  rotation  around  q,  and 
in  summary  multiplying  a  vector  x  on  the  right  by  a  quatern 
ion  q  (where  q  is  the  vector  part  of  q)  rotates  the  vector 
negatively  by  an  angle  tan”^-^-^  and  changes  its  length  by 
a  1 q ]  factor  (Ref  5). 
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1  rr}^ 

Now  consider  the  right  hand  quotient  of  q  and  x. 


q"1*  = 


qQx  +  (xXq)  +  (q*x) 
N(qj 


Since  x  and  q  are  perpendicular, 

q*x  =  0 

and  the  results  are  similar  to  the  previous  xq  except  the 
magnitude  is  reduced. 


“I”*! 

q  xl 


i7|Sl|g|\2  , 

(%wy] 

L\  NTqT/ 

\iTTq rv  J 

7  |  Xi  \2  !  |  I 

,11/2 

_\NlqT/  ^ 

J 

lxl |q| 

N(q) 

Since  q 


=  (q)^ 


1/2 


then 


N  (q) 


Therefore , 


q-1*! 


and 


q| 


tana  = 


!xl  I q | 

N(q) 

N(q) 


.*  q 


q0lxl 


=  JJlL 
qo 
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Because  x  is  perpendicular  to  q,  xq  is  also  perpendicular 
to  q  and  if  q™^  is  applied  to  xq  to  form  a  vector  y 

y  =  q  (xq) 

y  =  q”^xq  . 

The  vector  y  is  rotated  from  x  negatively  by  an  angle  of 
9  =  2a.  Now  consider  a  special  case  where  the  vector  x  lies 
along  q  and  has  a  magnitude  of  X[q| .  Then 

x  =  Xq 

where  X  is  a  real  number.  If  the  same  operation  is  carried 
out 

Y  =  q_1xq  =  q-1 (Xq) q 

y  =  q^txfq  -  qQ)]q 
=  X [q"1q  -  q_1qQJq 
=  X (1  -  qQq"1]q 
=  X[q  -  q0q_1q] 

=  X [q  -  qQJ 
=  Xq  =  x 

So,  a  vector  that  lies  along  q  and  operated  upon  by  q  1xq  is 
not  changed  in  magnitude  or  direction  (Ref  5) . 
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Finally ,  consider  a  general  vector  x  and  a  quaternion 
q  =  q^.  +  q.  The  vector  x  can  be  resolved  into  components 
(Ref  5) 


x  =  xx  +  x2 


along  q  and  perpendicular  to  q  such  that 


where 


^  _  q.x  -k 

1  |q|2  q 


and 


where 


q»x 

151 


x |  along  q 


x„  = 


(qXx)  X  q 

|q|2 


x|  perpendicular  to  q 


Fig.  5.  Components  of  x 

The  dot  product  of  q  and  x  (oriented  as  in  Figure  5) 
with  the  angle  B  between  them  is 


q*x  =  I q  |x|cos6 
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or 


cos  8  = 


q»* 


q|  I* 


The  cross  product  of  q  and  x  is 


or 


Therefore, 


qXx  =  j q| | x| sin8 


sin8  = 


»  _  qXx 


q|  ix 


X, 


x  cos8  “  x 


JH2L 


q| 


or 


and 


I  q| 


1*1  sine  =  |xj 

|q  lx 


or 


q 


Now  a  vector  y  exists  such  that 


-1-* 

y  =  q  xq 


—  1 

q  (x^  +  x2)q 

-1-k  -!-» 

q  xxq  +  q  x2q 

*1  +  'J'1x2q  =  Yl  +  Y2 
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where 


=  x^  along  q 

1  (  —A 

=  q  X2q  perpendicular  to  q 

The  vector  y2  is  normal  to  q  and  is  a  negative  rotation 
about  q  through  an  angle 

0  =  2  tan""^ 

q0 

When  q  is  oriented  along  the  Cartesian  x-axis,  the  rotation 
from  x2  to  y?  is  a  negative  rotation,  but  this  negative  space 
rotation  is  equivalent  to  a  positive  rotation  of  the  coordi¬ 
nate  system,  with  q^xq  representing  a  positive  rotation 
of  coordinates  and  qxq*'1  a  negative  rotation  (Ref  5) . 

Tilt  and  Roll  Quaternions 

To  utilize  quaternions  in  computer  programs  rotating 
coordinate  systems  in  three  dimensions ,  they  must  be  fac¬ 
tored  into  a  convenient  format  (Ref  15 .  This  format  factors 
the  quaternion  into  two  quaternions  which  divide  the  rota¬ 
tion  into  two  successive  transformations,  tilt  and  roll. 
Because  cf  this,  the  factored  quaternions  are  called  the 
tilt  quaternion  and  the  roll  quaternion  (Ref  5) . 

A  quaternion  q  =  +  q^i  +  q2j  +  q3k  can  be  written 

ss  q  =  qQ  +  q^i  +  k(q3  +  iq2)  ,  where  k(q3  +  iq2)  =  jq2  +  3cq3- 
For  any  complex  number  a  +  bi. 


and 

*2 
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k(a  +  bi) 


(a  -  bi)k 


9 


or  k  times  a  complex  number  is  equal  to  its  conjugate  times 
k  (Ref  1) .  Now  q  can  be  written 


q  =  Q 


01 


kQ 
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where 


®01  q0  +  iql  and  ^32  ~  q3  +  iq2  * 

The  conjugate  and  Norm  can  be  written: 

q*  =  °01  "  kQ32  '  where  °01  =  q0  ‘  iql 

N(q)  =  |Q01|2  +  Iq32|2  =  qo  +  <ii  +  +  q3 

Similarly,  a  vector  3c  «  ix1  +  jx2  +  kx3,  may  be  factored 
into 

x  =  ix1  +  kZ  where  Z  =  x3  +  ix2 

or 

x  =  (x.^  +  kZ)  i  where  Z  =  x2  -  ix3 

For  any  two  quaternions  q  and  r,  their  product  may  now 
be  written 

qr  =  (Qq1  +  kQ32) (RQ1  +  kR32) 

=  Q01R01  +  Q01kR32  +  kQ32R01  +  kQ32kR32 
=  Q01R01  +  kQ01R32  +  k®32R01  +  k  Q32R32 
=  (Q01R01  ~  Q32R32J  +  k(Q01R32  +  Q32R01) 
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Consider  a  quaternion  q  *  QQ1  +  kQ32  with  qQ  0.  The 
quaternion  can  be  factored  in  two  parts 

q  =  (1  +  k&Qoi 


=  (1  +  kT)QQ1  =  tr  where  T  =  ^ 


32 

01 


or  T  =  t3  +  it2 


The  tilt  quaternion  t  is  therefore 


t  =  1  +  kT 


and  the  roll  quaternion  is 


=  Q 


01 


=  q0  +  iqx 


Tilt 

The  vector  portion  of  the  tilt  quaternion  kT  is  located 
in  the  plane  normal  to  the  x-axis.  It  transforms  the  x-axis 
by  t-1xt  a  new  axis  x'  representing  a  tilt  from  x  results 


Fig.  6.  Tilt 
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The  vector  t  remains  normal  to  x  and  x '  and  represents  the 
axis  of  rotation  of  x  to  x '  through  an  angle  0 .  The  tangent 
of  0  being 


tan 


0 

1 


I  T  | 


♦  t§) 1/2  . 


The  y  and  s  axis  are  also  tilted  to  y 1  and  z '  axis  but  this 
is  not  shown  in  Figure  6  (Ref  5) . 


Roll 

The  roll  quaternion  r  has  the  form  r  =  (1  +  ir1> .  The 
inverse  is 


So  for  a  given  vector  x  = 
where  Z  =  x2  -  ix3  and  Z1 

x1 


Nlrf  =  * 

(x1  +  kZ)i  and  an  x’ 
*  x2  "  ix3 ' 

=  r_1xr  .  (Ref  5) 


(x£  +  kZ ! )  i 


r  represents  a  roll  about  the  x-axis  through  an  angle  <p . 
This  is  accomplished  by 

(x£  +  kZ' )i  =  r”1(x1  +  kZ)  ir 

or 


Thus 


(x£  +  kZ')i 


(x,  +  kZ  -)i 
1  r* 


x 


1 
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and 


*  Z 


r*r 


r*  rr 

The  angle  <j>  is  defined  by 


nTFT 


tan 


£  =  liL  = 


r  =  r. 


(Ref  5) 


Therefore 


r  =  1  +  i  tan  j 


rr*  =  1  +  tan 


2  4> 

I 


and 


d>  (h 

(1  +  i  tan  j)  (1  +  i  tan  j) 


*  1  -  tan  ( j)  +  i  2  tan  j 


cos<{>  +  i  sin<$> 


For  a  positive  roll 


Z  '  =  Z  (cos<l>  +  i  sin<j>) 


or 


x2  +  ix^  =  (x2  +  ix3)  (cos<j>  +  i  sin<M 

and 

x£  =  x2cos£{>  -  x3sin<J> 
x3  =  x2sim{>  +  x3cos<}> 


Therefore 


Z’  =  Ze1^  .  (Ref  5) 
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Appendix  B 

Convenient  Computer  Coding  Relations 


Complex  Factoring 

The  tilt  quaternion  t  represents  a  transfer  of  coordi¬ 
nates  from  an  unprimed  system  to  a  primed  system  (Ref  5) . 

A  general  vector  x  defined: 

x  =  (x1  +  kZ)i  with  Z  =  x2  -  ix3 

will  become 

x'  =  (x£  +  kZ 1 )  i  with  Z '  =  x2  ~  ^"x3 

This  is  accomplished  by: 

x*  =  t"1  x  t 

where 

t  =  1  +  kT  and  T  =  t3  +  it2 
By  expanding 

*•  ■  ?  (1  +  kT) 

=  (---  kTj,  x  (1  +  kT) 

i  +  It!2 

Let 


d 


2 

1  +  |  T  i  2 


(3) 
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then 


xs  *  |  (1  -  kT)  x  (1  +  kT) 

=  |  (1  -  kT)  (k±  +  kZ)  i  <1  *  kT) 

=  |  (1  -  kT) ( xl  +  kZ) (i  -  kT) i 


Expanding 

(xj  +  kZMi  =  §  [  (xx  +  kZ  -  kTxL  -  kTkZ)  (1  -  kT)  i] 

=  ^  [x^  -  xi5cT  +  kZ  -  kZkT  -  kTx^  +  kTx^kT  - 

kTkZ  +  kTkZkT] i 

=  2  fx^  “  kx^T  +  kZ  +  ZT  -  kTx^  -  x^TT  +  TZ  -  kT^Z]! 

=  |  [x1  +  ZT  -  xxTT  +  TZ  +  k(Z  -  2Tx1  -  T2Z) ] i 

(x£  4-  kZ ' )  i  -  |  [Xi  -  x1lTj2  +  ZT  +  TZ  +  k(Z  -  2Tx1 

-  T2Z)  ]  i 

x£  +  kzs  =|  [Xx  -  xjTl2  +  ZT  +  TZ]  +|  k [Z  -  2Txl 
-  T2Z] 


Therefore 


I  tx1(l 


T| 2)+  ZT  +  TZ] 


and 


=  f  [Z 


2Txx  -  T2Z] 


(Ref  5)  (4) 


Let 


b 


x. 


+  x£ 
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then 


b  S  X1  +  f  txl^  ~  iT!2)  +  2T  +  TZ] 

but 

ZT  =  (x2  +  ix3)  (t^  +  it2)  =  x2b3  “  x3b2  +  ^■^x2t2  +  x3b3^ 

and 

TZ  =  (t^  -  it2)  (x2  -  ix^)  =  x2t3  "  x3t2  ~  i(x3t3  +  x2fc2^  • 

So 


ZT  +  TZ  =  2x2t3  -  2x3t2 

and 

b  =  |  [|  xx  +  x1U  -  | T |  2)  +  2Re  (TZ)  3 

where 

Re(TS)  =  x2 t3  -  x3t2 

Substituting 


b  =  |  [x^l  +  |  T  |  2 )  +  Xx(l  -  |  T  |  2)  +  2Re  (TZ)  3 

b  =  |  [2xx  +  2Re(TZ)3 


=  d  [x-j^  +  Re(TZ)  ] 

(5) 

=  3  [xx  +  t3x2  -  t2x3]  . 

(Ref  5) 

(6) 
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From  (4)  we  have 


Z'  =  |  [Z  -  2Txx  -  T2Z] 


(7) 


Subtracting  Z  from  both  sides 


Z '  -  Z  = 


|  [Z  -  2Txl  -  T2Z]  -  Z 


=  f  Z  -  Z  -  f  [2Tx.  +  T2Z] 


=  Z(|  -  1)  -  |  [2Txx  +  T2Z] 


but 


_  -  1  = 


l  +  It 


-  t 


1  +  T 


I  (-TT) 


Substituting 

Z'  -  Z 


|  ZTT  -  |  [2Tx1  +  T2Z] 


=  -  |  T[2xx  +  ZT  +  TZ] 
=  -  |  T[2x;l  +  2Re (TZ)  ] 


Z  -  Z'  =  - 


bT 


So 


Z 1  =  Z  -  bT 


(Ref  5) 


(8) 


A  reversal  of  the  tilt  transformation 


x'  =  t”1  x  t 


yields 


t  x1  t"1  =  tt-1  x  tt"1  =  x 
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This  leads  to  the  equations  of  the  inverse  transformation 


d 


2 

l  +  |t|  2 


(9) 


b  =  d(x£  -  *  fc2x3^ 


(10) 


x.^  =  b  -  x^  (11) 

Z  =  Z'  +  bT  .  (Ref  5)  (12) 

Forming  the  Tilt  for  Range 

For  a  given  vector  r  in  an  unprimed  coordinate  system 
there  exists  a  quaternion  t  which  will  tilt  the  vector  r 
to  the  x'-axis  (Ref  5).  Let  r  be  defined  by: 

r  =  (r,  +  kZ) i  where  Z  =  r~  -  ir_ 

j.  2  3 

and  |r|  =  R.  Now  to  determine  t  which  will  tilt  the  x'-axis 
to  r  by 

r'  =  t  r  t 


or 


(r£  +  kZ')i  =  t"1  r  t  . 

But  rj^  »  R  and  Z'  =  0  since  r  is  to  be  oriented  along  the 
x'-axis.  Thus  utilizing 


and 


b  *  r^  +  r£  =  r1  +  R 

Z'  =  Z  -  bT  =  0 
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v-  ,  > 


the  tilt  quaternion  t  =  1  +  kT  can  be  found* 


bT  *  Z 


z  _  r2  “  ir3 

T  “  b  ■  -mr 

If  on  the  other  hand,  a  quaternion  t  is  desired  which  will 
tilt  the  x'-axis  to  r,  the  following  equations  are  used 
( |  r |  =  R,  Ze  =0,  r =  R)  . 


b  =  d[r£  -  Re(TZ') ]  =  dR 


where 


where 


d  = 


1  +  T 


rl  *  b  -  R  =  R(d  -  1) 


Z  =  bT  =  dRT 


d  -  1  = 


d  -  1  = 


1  +  |  T I 
1  -  tan2  | 
1  +  tan2  | 


=  J 


cos2  (^-) 


d  -  1  =  cos0 


=  R (d  -  1)  =  R  cos0 
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Roll  Representation  Due  to  Tilt 

Because  missile  seeker  heads  are  gimbaled  to  move  in 
azimuth  and  elevation  Euler  angles,  a  tilt  of  the  x-axis 
will  be  accomplished  by  movements  constrained  by  the  fixed 
axes  (Ref  5) .  The  tilt  quaternion 


t  =  1  +  kT 


where 


T 


r2  *  ir 
R  +  r1 


3 


+  ifc2 


transforms  from  i,  j,  k  to  i',  j',  k'  coordinates  by  tilting 
i  through  an  angle  t .  The  case  illustrated  in  Figure  7 
shows  a  change  in  azimuth  followed  by  a  change  in  elevation 
(Ref  5)  .  To  model  this  motion,  the  desired  roil  angle  <j>  is 
assumed  known  and  through  the  inverse  of  this  roll  transform 
the  coordinates  to  the  tilt  position  defined  by  t-  Then  by 
using  t  1  the  original  position  is  obtained  (Ref  5) . 
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The  j  coordinate  in  the  final  position  following  the 
azimuth  and  elevation  change  can  be  defined  in  factored  form 


j 


s 


+  k  js)i 


where 


js  =  01  +  js2^  +  0k  • 

Since  jg^  =  0  and  ki  =  j 


J 

s 


1 


For  an  inverse  or  negative  roll 


Z' 


or 


j,  =  j  e-i4>  = 

s  s 

3si  "  iii  -  0 

and  for  an  inverse  tilt 

t  =  (1  +  kT) 

T  =  (t3  +  it2) 


where 


d 


2 

1  +  |  T  j  2 


(Ref  5) 
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and 


b  =  d [ jgX  ”  Re(TJg)  ]  =  -  dRe(Te"l4)) 

Therefore 

J"  =  J'  +  bT  =  e-1^  -  dRe(Te“ic^)T 

b  b 

bat 


Js  =  (js2  ’  i^s3)  where  ^ "3  =  0 

So 


-  ImJ"  =  0 

s 

and 

Im[e-1^  -  dRe (Te-1^) ] T  =  0  . 

Since  for  any  complex  number  A 

A--~-A  =  Im(A)  and  A  |-A  =  Re (A) 

~  [{e"1(|)  -  dRe  (Te-1<^)  T)  -  (e"^  -  dRe (Te-1^) T) 

Collecting  terms 

e"1^  -  dRe (Te""1^)  T  =  eitf)  -  dRe  (Te“1(j!)  S 
e-1^  -  e3^  =  dRe(Te''i<^)  (T  -  T) 

e"1^  -  e1^  =  2idRe (Te”1^) Im{T) 

e“1<j>  -  e1(J)  =  idlm(T)  (Te-1^  +  Te1^) 


.  (Ref  5) 


(Ref  5) 
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Multiplying  both  sides  by  e1^ 

1  -  e21^  =  idlm(T)  (T  +  Te2l<h 

Collecting  terms, 

1  -  idlm (T)  T  =  e2l<{>[l  +  idIm(T)T}  .  (Kef  5) 


Let 


A  =  1  -  id!m(T)T 


then 


A  =  1  +  idlm(T) T 


and 


_ 2icf)  _  1  -  idIm(T)T  _  A 

1  +  idIm(T)T  A 


e2i* 


A 

A 


(Ref  5) 


Since 


r  = 


2  1/2 


I  A  i 


(13) 


The  roll  quaternion  r  acts  through  the  angle  4>  to  restore 
the  tilted  coordinates  to  the  zero  roll  position  where  the 
tilt  quaternion  transforms  them  to  their  original  position. 
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a  Roll  for  Range 
The  roll  transfer  relation  is  formed  from  the  relations 


r 

r* 


A 

I A I 


(14) 


This  is  accomplished  by  utilizing  the  relationships  devel¬ 
oped  previously  for  a  tilt  for  range  where 


t  =  1  +  kT  T  = 


r2  “  ir3 

R  +  r. 


and 


b  =  dR 
y_  R  +  r. 

d  =  ^  i 

R  R 


Then 


A  =  1  -  idlm(T) T  =  1  -  i 


+  rl  /  ~r3  \/r 2  +  ir3\ 


.2  2 

-r-r,,  -  xr,  r^  -  r^r^x 

A  =  1  -  i  ■  -i--:— r  =  1  - 


R(R  +  rx) 


and  finally 


A  =  1  - 


R  (R  + 


r2r3 


R  (R  +  rx)  +  R(R  +  r1)  1 


Rotating  Coordinate  Systems 

Let  two  rotating  coordinate  systems  be  related  by  a 

tilt 


t  =  1  +  kT 


that  is 
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X'  =  t  X  t 


with  angular  velocities 


<o  =  ico^  +  jWj  +  kco. 


and 


a) 1 


iw^  +  jco^  +  kt0^ 


at  a  later  time  tQ  +  At  the  systems  will  relate  through  an 
angle  jwjAt  about  the  w  axis.  Therefore,  the  cu  rotation 
is  similar  to  a  quaternion  (Ref  5) 


Sit  =  1  +  “q 


or 


,  CO  I  CO  I  At 

q, .  =  1  +  -=—  tan  i  !-■■■■ 

At  i  —•]  2 

co 


Since  for  small  angles  tana  =  a 


Sit  =  1  +  T4  (Ref  5) 


coAt  ,  -*  ,  j 

Let  — 2~  equal  a  vector  x,  and  x 

of  x  can  be  defined 


ix^  +  kX. 


The  components 


co-^At 


<02it 


and  x. 


co^At 


and  let 


X  = 


(CO2  +  ico2)At 


OAt 

2 


0  =  co^  +  icoj  (Ref  5i 
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There  also  exists  a  q^t  such  that 


q£t  -  1  +  ix^  +  kX*  =  1  +  . 

Let  s  be  a  quaternion  that  relates  the  unprimed  coordinate 
system  to  the  primed  system  at  time  (tQ  +  At) .  Assuming  s 
is  a  tilt  as  was  t  then 

s  *  1  +  kS  =  1  +  k(T  +  AT) 

Now  the  quaternions  qAts  and  tq^t  both  transform  from 
unpriraed  to  primed  coordinates  at  time  t  (Ref  5) .  They 
also  transform  from  unprimed  to  primed  coordinates  at  a 
t  +  At  and  are  therefore  equivalent  (Ref  5) 

qAts  =  atq£t 

where  a  is  a  real  number.  Let  a  =  1  +  h  so  that  when 

At  =  0,  h  =  0  and  q^fcs  =  tq^fc.  Substituting  for  qAfc/  s,  t 

and  q^fc  and  neglecting  second  order  terms 

LHS  =  (1  +  ix1  -  XT)  +  k(X  +  T  +  AT  -  ix-jT)  +  ord(At2) 

RHS  =  (1  +  ix£  -  TX'  +  h)  +  k(X*  +  T  +  ix£T  +  hT)  +  ord(At)2. 

o  2 

where  ord(At)  is  on  the  order  of  At  .  Collecting  corres¬ 
ponding  parts: 


ix1  -  XT 

*  ix|  -  TXf  +  h 

(15) 

X  +  AT  -  ix^T 

*  X*  +  hT  +  ixjT 

.  (Ref  5)  (16) 

Equating  the  real  and  imaginary  parts  of  (15) : 

h  =  Re  (TX '  -  XT)  (17) 
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x^  =*  Im(XT)  **  x^  +  Im(TX) 

*  xX  “  Im(TX' )  (Ref  5)  (18) 

If  a  quaternion  p  is  now  defined  as 

p  =  (x£  -  Ira(TX' )  ~ 

then  p  =  tu£  -  Im(TO')  =  oj^  +  im(TO)  which  constrains  w  and 
w'  so  that  p  will  remain  a  tilt  from  unprimed  to  primed 
coordinates  (Ref  5) .  From  (17) 

h  =  Re  (TX ’  -  XT) 

=  Re (TX 1 )  -  Re (XT)  =  Re(TX')  -  Re (TX) 

thus 

h  =  Re[T(X’  -  X)]  .  (Ref  5) 

Mow  from  (16) 

AT  =  X'  -  X  +  T[h  +  Kac^  +  x£)  3 

*  (X‘  -  X)  +  T  [Re  (T  (X '  -  X))  +  i(x1  +  xp] 

So 

AT  *  (1  +  j  T | 2)  (X*  -  X)  +  T[Re(T(X'  -  X))  -  T(X’  -  X) 

+  i  (x^  +  x£)  3 

«  \  (X'  -  X)  +  iT[xi  -  Im (TX ' )  +  X,  +  Im(TX) ] 

a  11 
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and 


AT  ^  |  [X'  -  X  +  2iT(x£  -  Irafx')]  (Ref  5} 
Dividing  both  sides  by  At  and  taking  the  limit  as  At  -*•  0 

T  *  j  (O'  -  O)  +  iTp  (Ref  5) 

where 

dx2  dx« 

O  =  <u3  +  i«2)  “2  ’  -3T  '  “3  “  TF 

As  long  as  p  =  +  Im  (TO)  =  tu^  -  Ira  (TO1)  the  unprimed  and 

primed  coordinates  will  remain  related  by  a  tilt.  (Ref  5) 

1  +  kT 

and  T  can  be  determined  by 

T  =  /  T  dt  .  (Ref  5) 

If  one  of  the  coordinate  systems  is  fixed  then  tu  =  0,  O  =  0 
and  p  =  0  giving  =  Im(TO')  and 

T  =  ~  .  (Ref  5)  (13) 

In  both  cases  T  is  numerically  integrated  to  find  the  tilt 
t  (Ref  5) .  To  relate  a  fixed  axis  system  to  a  non-rolling 
system  with  rolling  components  in  the  j'  k'  plane  only,  an 
intermediate  rolling  coordinate  system  R  with  angular  velocity 


ua  +  kO ' 
R1 


with  to-,  =  Im(TO')  so  that 
R1 
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T 


O' 

T  * 


Therefore,  the  tilt  relating  R  to  the  fixed  (F)  system  can 
be  found  by  numerical  integration.  If  the  quaternion  q(t) 
relating  R  to  F  is 

q  ( t)  =  1  +  kT(t) 

where  T(t)  =  T  +  f  T  at  and  a  roll  quaternion  r(t)  relating 
the  non-rolling  system  to  the  rolling  system: 

r (t)  =  (fc)/2  , 

the  quaternion  s(t)  relating  the  non-rolling  to  the  fixed 
system  is: 

s  { t)  =  q { t)  r  (t)  =  [1  +  kT(t)]e1<|)(t)/2  .  (Ref  5) 

Therefore,  a  <Mt)  must  be  found  such  that  of  the  non¬ 
rolling  system  equals  zero  (Ref  5) . 

The  quaternion  transforms  s(t)  from  F  to  non-rolling 

(NR)  axes  at  time  t.  For  a  small  time  interval  At  and  a 

-1 

time  (t  -  At),  the  quaternion  [s(t  -  At)]  transforms  NR 
to  F  and  if 

p  *  [s(t  -  At)]"1  s { t)  =  PQ1  +  kP32 
transforms  NR  at  (t  -  At)  to  NR  at  time  (t) .  Since  uiNR  =  0 

<u>NR^)  (At)  =  Im(P01)  =  0 

substituting  into  p 
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~  At)/23[1  ,  kT(t  _  At)  ]  [1  +  kT(t)  3ei<i)(t)/2 

Njg(t  _  At)] 


and 


P 


01 


[e  lsi(t  At)/2]  [i  +  :r  ( t  ~  At)  T(t)  le1^  (t)/2 

— — - - - - -  7""~~ -  ~  2  - — ~  -  (Ref  5) 

1  +  T(t  -  At)  * 


Let 


A<f>  =  <f>(t)  -  <j>  (t  -  At) 

then 

P01  =  El  +  T(t  -  At)  T(t)  ]elA4,/2  .  (Ref  5) 

If  Im{P_ , )  =  0  then 
0 1 

arg  [1  +  T(t  -  At)  T(t)  +  arg  (elA(i)/2)  ]  =  0  CRef  5) 

If  A  is  a  complex  number,  then 

A  =  a  +  bi 


When  expressed  in  polar  form 


U2  H-b2)1/2 


LCa2  +  b2)1''2 


+  i 


1(aTTW7I. 


or 


A  =  rCcoscJ)  +  i  sin<|))  =  re^ 

where  r  *  Ca2  +  b2)^2  and  cos<$>  =  sin<J>  =  p  (Ref  9). 
The  argument  (arg)  of  A  is: 

argA  »  $  .  CRef  9) 
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Therefore  if  b(Im&)  is  equal  to  zero  then  cos#  *  1  and 
sin#  *  0.  Which  means  that  #  ®  0. 

Since  consists  of  two  complex  factors 

[1  +  T(t  -  At')  T  (t)  ]  and  e 

their  arguments  must  sum  to  zero,  otherwise  #  ^  0.  Then 
arg[l  +  T(t  -  At)  T(t)J  =  -arg telA4>/2) 

but 

arg  ei4^2  =  ^ 

and 

=  -arg[l  +  T  (t  -  At)  T(t)J 

A#  *  -2arg[l  +  T  (t  -  At)  TCt)J 

or  A#  *  -2arg[l  +  T(t  -  At)  TTtf]  C20) 

The  model  numerically  integrates  (.20)  using  the  present  and 
past  values  of  T  which  is  numerically  integrated  from  equa¬ 
tion  (19)  (Ref  5)  . 


Consider  a  coordinate  system  as  in  Figure  8.  The  vec¬ 
tor  kT,  where  T  -  t^  +  it2  in  the  complex  Xg/X-,  plane, 
defines  a  new  coordinate  system  which  is  rotated  about  the 
x,  axis  through  an  angle  y  with  major  axis  y~  ,  y^,  y3- 


Fig.  8.  x^  Rotation 

Now  consider  a  rotation  about  the  y3  axis  as  depicted  in 
Figure  3.  These  two  rotations  represent  a  tilt  of  the  x^ 
axis  to  the  y £  position  in  space. 


To  accomplish  the  tilt  from  x^  to  x’  through  the  angle  0, 
using  direction  cosines,  a  two  step  process  must  be  utilized. 
Let  the  rotation  about  x^  be  represented  by  direction  cosine 
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matrix  A  and  the  rotation  about  y^  be  represented  by  direc¬ 
tion  cosine  matrix  B  (Ref  4) .  Then 

y  -  Ax 

rotates  x  to  y  in  x  coordinates , 

Y*  =  By 

rotates  y  to  y'  in  y  coordinates.  So  the  original  vector  x 
rotated  by  A  then  B  is 

x '  =  ATBy 

in  x  coordinates.  (For  an  orthogonal  matrix  A  =  A  )  The 
completed  tilt  can  be  represented  by 

x'  =  Cx 


where 


C  =  A  BA  . 


From  Figure  8 


siny  =  -  — r—  ,  and  cosy 

It| 


From  the  quaternion  definition  of  t. 


tan  2  = 


then 


x  =  2  tan’1|T| 
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1 


Fig. 10.  tan~1|T| 


From  Figure  10 


cos (tan”1! T|) 
sin (tan”1 | T| ) 


(1  +  |T|2)1/2 

IT] 

77T  I  m  I  2  \  1/2 


COST 

COST 


and 


sinT 


cos  2 (tan”1 [ T | ) 

2  cos2 (tan”1! T| )  -  1 

2  _  it  1t1 2  _  l  -  It] 2 

1-4-  1 T I  2  1  +  [Tj  2  1  +  |T]  2 


sin  2 (tan”1 | T | ) 

2  sin (tan  ^|t|)  cos (tan  1 ( T | ) 


I  T  { 

1 

_  (1  +  1 T  j  2)  1//2_ 

.  (1  +  |  T  |  2)  1/2_ 

2  T| 

—  2 
1  +  T  c 


Now  that  the  functions  siny*  cosy,  sinT  and  cost  are  defined, 
the  transformation  matrices  can  be  formed: 
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1 


0 


0 


and 


A 


0  cosy  siny 
0  -siny  cosy 


B 


cost  sinx  0 

sinx  cost  0 

0  0  1 


c  =  atba 


c  = 


COST 

sinx cosy 
siny sinx 


cosysinx 

2  2 
cos  ycosx+sin  y 

sinycosycosx-sinycosy 


siny sinx 

sinycosycosx-sinycosy 

2  2 
sin  ycosx+cos  y 


Substituting  and  carrying  out  the  vector  multiplication 

1  —  ] <r|  2  ^t,  ^t- 

x«  =  ± - LLL_  x,  +  - x~  +  - x, 

1  1  +  It]2  1  l  +  1 t | 2  2  i  +  It]2  3 


2t 


k2  =  " 


3  1  +  |T|2  -  2t2  2t,t, 

X,  +  - s - -  X~  +  - — — —S'  X, 

*  1  .  !  m  1  <£ 


1  +  I  T I  2  1  1  +  I T  ]  2 


1  +  {  T  j 


2t 


X3 


,  2t~t,  1  +  | T I  2  -  2t2 

x,  +  - — — — s-  x„  +  - 2 — “  x3  • 


l  +  It!2  1  l  +  ] t 1 2  2 


i  +  t! 


From  the  tilt  derivation  for  a  positive  rotation  of  coordi¬ 
nates 


xj^  =  b  -  x^  and 


Z’  =  Z  -  bT 
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or 


xJ,  -  ix-j  -  (x2  -  ix3)  -  b(t3  +  it2) 

therefore 


and 


where 


x3  =  x3  +  bt2 


b  =  d(x,  +  t,x~  -  t-x.,)  and  d  =  - 5- 

13223  1  +  |T|Z 


Substituting 


x 


x2 


xl  = 


_  1  ~  !  T  | 


2t- 


2t. 


x,  + 


x„  - 


1  +  | t | 2  1  1  +  It!  2  2  1  -h  It!2  3 


2t. 


1  +  |T 
2t„ 


1  +  |T|2  -  2t2  2t2t3 

21  1  +  |  T  |  2  2  1  +  I T 1 2  3 


1  +  I T I  2  1  1  +  !  T } 


x,  + 


2t2t3  1  +  | T | 2  -  2t2 

2  X2  +  ,  ,  ,„,2  X3 


1  +  T 


Therefore,  the  tilt  quaternion  gives  the  same  result  as  the 
three  matrix  multiplications . 

Similarly  for  a  negative  roll  but  a  positive  coordinate 


transformation 


xs  =  Ax  where  A 


10  0 
0  cos0  sin0 
0  -sin0  cos0 
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and 


x£  =  x2cos<?>  +  x^simj) 
x^  =  XpSin<|)  +  x^costi? 


From  roll  quaternion  development 


r  =  1  -  i  tan  j 


2  0 

1  +  tan^  j 


r 

.2 


2  cf? 

1  -  tan^  j  ~  i  2  tan 


=  cos<fc  -  i  sin4> 


$ 

1 


Then 


z*  =  (cos<|>  -  i  sin<|>)  (x2  +  ix3) 


x2  +  ix^  =  x2cos4>  +  x3sin®  -  ix2sin$  +  ix3cos0 


or 


x2  =  x^ostf1  +  x3sin6 
x^  =  -  x2sin<l>  +  x3cos® 
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Coordinate  Systems 


From  the  discussion  on  rolling  coordinate  systems  in 
Appendix  B,  a  p  quaternion  is  formed  by: 


P  =  qAt  3 

where  p  relates  a  rolling  and  non-rolling  set  of  coordinate 
systems  at  any  time  t.  From  the  development  of  p,  equation 
C19)  is  derived  which  implies  the  parameter  of  the  rota¬ 
tional  velocity  vector  w  is 

u)£  =  Im  CTO  * ) 

This  is  equivalent  to: 


=  ImCt3  -  it2)  Cw*  +  iw2) 
=  t3u)£  -  t2w3 


C21) 


To  arrive  at  C21)  utilizing  conventional  direction 
cosines,  it  must  first  be  understood  that  the  product  of  two 
quaternions  will  result  in  a  negative  coordinate  rotation. 
This  is  a  reversal  of  the  direction  discussed  previously  in 
this  appendix  and  result  in  the  following  changes: 


siny 


,  cosy 


COST 


1  +  |T|2 


sinx 


1  +  |T|2 


In  addition,  the  cosine  matrices  A  and  B  are  now 
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and 


1  0 


0 


( 


A  " 


0 

0 


cosy  -sinY 
sinY  cosy 


B 


cost  -sinx  0 

sinx  cost  0 

0  0  1 


To  describe  the  vector  to  which  relates  the  non- 
coordinate  system  to  the  rolling  coordinate  system, 
following  relation  is  used: 

wdCt)  =  CCt) 


where 


C  =  ABA 


The  vector  53  can  be  found: 


where 


or 


o'  =  CCt)CTCt) 


CCt)  «  ^  (ABAT) 


.m  .mm  .m  m 

w  ®  AA  +  ABB  Aa  +  ABCA  A)  (AB)  i 


Substituting  for  A,  A,  AT,  AT,  B,  B  and  Bl 


rolling 

the 
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0 


YsinysinT-TcosY  -YCOsysinT-Tsiny 


to 


-YsinysinT+Tcosy  0 

YCOSYsinx+xsinY  y-ycosT 


-y+YCOST 

0 


which  implies 


-i  - 

y(l  -  cost) 

“2  - 

•  • 

-  ycosysinx  -  xsiny 

and 

“3  - 

-  ysinysinx  +  xcosy 

/  /  • 

Utilizing  the  equations  for  an^  ^3*  Y  can  be  solved  for: 

.  to^cosy  +  oo^siny 

Y  -  ”  sint  * 

Substituting  into  the  equation  for 


0) 


^  =  (w^cosy  +  co^siny)  C 


cost  -  1, 
sinx 


Substituting  for  cosy, 


(.  ^ 

1  V  |T| 


siny,  cost  and  sinx: 

^2^  /  2  f  T{  2  \  fl  +  [T{ 

T|  /  \  1  +  1 T (  /  \  2  ( T I 


or 

“l  “  “2fc3  "  u3t2  (22) 

Equations  C21)  and  (22)  are  identical  and  the  same  results 
are  obtained. 
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Appendix  D 

Computer  Listing  and  Guide 


introduction 

The  subroutine  SIGNAL  was  modified  so  that  a  table  of 
signal  strengths  for  various  aspects  and  ranges  could  serve 
as  a  data  base  for  calculations  of  maximum  seeker  line  of 
sight  (LOS)  rates  for  target  tracking.  To  accomplish  this, 
azimuth  and  elevation  data  arrays  (Figure  11)  were  inserted 
along  with  the  proper  code  (Figure  12 )  to  evaluate  the 
geometry  and  extract  the  proper  values  from  the  data  base 
to  interpolate  a  proper  signal  strength. 

The  missile  model  uses  many  common  variables  so  that 
quaternions  need  to  be  defined  once  and  then  through  the 
common  variables  be  available  to  various  subroutines  as 
necessary.  The  program  DRIVER  manages  the  overall  model  by 
reading  the  input  data,  initializing  the  variables  and  set¬ 
ting  up  quaternion  relationships  among  the  various  coordi¬ 
nate  systems.  Where  appropriate,  the  coding  from  DRIVER 
has  been , extracted  to  help  define  the  work  variables  which 
form  the  tilt  and  roll  quaternions  utilized  in  SIGNAL. 

The  following  sections  describe  the  code  used  to  modify 
SIGNAL.  The  reader  is  encouraged  to  read  Appendices  A  and  B 
prior  to  reading  the  Azimuth  and  Elevation  section.  The 
final  section  is  an  example  problem  to  demonstrate  how  the 
interpolation  technique  works. 
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Pig.  12b.  Code 


Range 


The  first  computation 

ITEMPA  =  R  *  0.001  +  1 

takes  the  actual  magnitude  of  the  range  vector  and  divides 
it  by  1000.  This  indexes  the  range  into  a  range  bin  with 
respect  to  the  data.  The  data  arrays  are  constructed  such 
that  the  row  value  corresponds  to  the  range  bin  and  the 
column  value  corresponds  to  the  LOS  azimuth/elevation  meas¬ 
ured  from  the  target  back  to  the  missile  seeker.  There  are 
four  range  indexes  and  13  angle  indexes. 

The  value  IRNG  is  next  assigned  by  determining  if 
ITEMPA  is  between  zero  or  one  (IRNG  =  1)  ,  between  two  and 
six  (IRNG  *  2)  ,  or  greater  than  six  (IRNG  =  3) .  These  val¬ 
ues  correspond  to  ranges 


IRNG  = 

1 

0  £ 

range  < 

1000 

IRNG  = 

2 

1000 

<  range 

<  5000 

IRNG  = 

3 

5000 

a> 

tr* 

9 

M 

VI 

<  13000 

Azimuth  and  Elevation 

The  first  step  in  setting  up  the  geometry  of  the 
missile- target  relationship  is  to  transform  the  line  of 
sight  (LOS)  from  the  fixed  (launch  aircraft)  coordinates  to 
U  (range)  coordinates.  This  is  done  by  calling  the  sub¬ 
routine  RESOLV  which  utilizes  a  direction  cosine  matrix 
(FTOU) .  The  calling  sequence  is  CALL  RESOLV (DLRF1,  DLRFR, 
-DLRFI ,  FTOU,  TEMP A,  ZWORKA) .  DLRF1  is  the  x  component. 
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DLRFR  is  the  y  component  and  DLRFI  is  the  z  component  of 
the  range  vector  in  fixed  (F)  coordinates.  (DLRFI  is  negative 
to  aid  later  computations.)  FTOU  specifies  the  direction 
cosine  matrix  to  be  used  (F  coordinates  to  U  coordinates) . 
Finally  TEMPA  and  ZWORKA  are  dummy  arguments  which  return 
the  range  vector  in  U  coordinates  in  factored  format 

r  =  (r1  +  kZ)i 

where  r^  is  the  TEMPA  value  returned  and  Z  =  ZWORKA  = 

{r 2  -  ir3).  The  components  r2  and  r 3  are  the  y  and  z  values 
of  r  in  the  U  system. 

Now  that  the  LOS  is  defined  in  the  U  system  which  con¬ 
tains  both  the  missile  and  target,  the  LOS  is  redefined  in 
the  target  system  to  determine  the  azimuth  and  elevation  of 
the  seeker  from  the  target  reference. 

Because  the  seeker  is  a  gimbaled  system,  both  a  tilt 
and  roll  are  required  to  match  the  LOS  from  the  target  to 
the  seeker  reference  system.  The  required  transformation 
utilizes  the  tilt  quaternion  t  =  1  +  k(ZTUT)  and  the  roll 
transformation  ZRUT. 

The  tilt  quaternion  is  constructed  first  through  the 
work  variables  TEMPA,  TEMPB,  WORKAR,  and  WORKAI: 

TEMPB  =  DUT* (TEMPA  +  WORKAR  *  TUTR  +  WORKAI  *  TUTI)  . 

DUT  is  defined  in  the  DRIVER  as  the  factor  d  of  the  tilt 
quaternion  (ZTUT)  from  U  to  T  coordinates  where 

DOT  -  - - - =r  .  (d  from  U  to  T) 

1  +  |T| 2 
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TUTR  is  the  real  part  and  TUTI  is  the  imaginary  part  of 
2TUT  which  is  also  computed  in  DRIVER.  The  sequence  from 
DRIVER  is  as  follows: 


DUT 


TEMPA 
TUTR 
TUTI 
TEMPA 
=  TWO/ 


=  VT  +  VTUl 
=  VTUR/TEMPA 
=  VTUI/TEMPA 
-  TUTI**2 

(ONE  +  TUTR* *2  +  TEMPA) 


The  variables  used  here  are  VT ,  the  magnitude  of  target 
velocity;  VTUl,  the  x-coordinate  of  velocity;  VTUR,  the 
y-coordinate  of  velocity;  and  VTUl,  the  z-coorainate  of 
velocity.  The  target  coordinate  system's  orientation  is 
defined  by  VTUl.  TEMPA  in  this  case  is 


TEMPA 

> 

II 

+ 

V1 

VTUl 

TUTR 

V2 

VTUR 

VT 

+ 

TUTI 

v3 

VTUl 

VT 

+ 

vT 

TEMPA  is  then  redefined 

V3 

TEMPA  =  - - - 

(VT  +  VL) 


Since  the  tilt 


t  =  (1  +  kT) 

V3 

and  T  =  t3  +  it2  where  ^  anc*  tg 

then  the  magnitude  of  T: 


VT  +  v. 
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DUT  therefore  is 


2  2 
1  +  +  tj 


1  •+  T 


Returning  to  the  SIGNAL  sequence,  TEMPB  is: 


TEMPB  =  b  =  d[r^  +  t3r2  -  t2r3]  , 


and 


rx  =  TEMPA 
r2  =  WORXAR 

-r3  =  WORKAI  (From  the  RESOLV  calling  parameter) 

t3  =  TUTR 
t2  =  TUTI  . 

With  these  parameters  the  new  range  vector  r '  can  be  con 
structed. 

* 

TEMPO  =  TEMPB  -  TEMPA 


is  equivalent  to 


b  -  r^  (from  App.  B  (11)) 


WORKER  ^  WORKAR  -  TUTR* TEMPS 


r2  t3b 


and 


WORKAI  =  WORKAI  -  TUTI*TEMPB 


=  -r,  -  t,b 


ZWORKB  is  defined  by 


Z WORKS  =  WORKAR  +  i WORKS I 


or  equivalently 


Z'  =  r2  “  ir2  •  (from  App.  : 


From  the  discussion  of  a  quaternion  times  a  vector 


Z1  =  z  -  bT  , 


where 


Z  =  r2  -  ir3  =  WORKAR  +  i (WORKAI)  . 


This  can  further  be  broken  down  into 


Z'  =  r,  -  bt^  -  i(r,  +  bt0) 


or 


WORKER  =  r.,  -  bt  , 


works  I  =  -r,  -  bt-, 


(8)) 
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The  roll  due  to  tilt  is  next  considered  using  the  roll  trans¬ 
formation  from  U  to  T  coordinates  (ZRUT)  which  is  defined  in 
DRIVER  by  the  sequence: 


WORKAR 

=  ONE  -  DUT*TEMPA 

WORKAI 

=  -DUT*TUTI*TUTR 

TEMPA 

=  CABS  ( ZWORKA) 

RUTR 

=  -WORKAR/TEMPA 

RUTI 

=  -WORKAI /TEMPA 

This  sequence  is  identical  to  the  development  of  the  roll 
due  to  tilt  discussion  where  the  roll  quaternion  r  =  (1  +  ir^) 
and  the  roll  transformation  is 

_  JL  —  _A_ 
r*  I A  |  • 

The  complex  number  is: 

A  =  1  -  idIm(T)T 

where  T  =  t,  +  i t_,  T  =  t_,  -  it-  and  d  =  DUT.  ZRUT  which 
is  equivalent  to  ~  is  found  by  — "  '  * 

WORKAR  =  Re  (A)  =  1  -  dt^ 

WORKAI  =  Im (A)  =  -dt2t3 

A  =  ZWORKA  =  WORKAR  +  i WORKAI 

therefore 

A  =  1  -  d(t2  +  it2t3^ 
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where 


iIm(T)T  = 


-  it0) 


+  it2t3 


Since 


and 


TEMPA  =  CABS  ( ZWORKA)  =  Ja  , 


RUTR  =  - 


Re  (A) 


RUT I  =  - 


Im(A) 


then  ZRUT  = 


The  roll  takes  place  about  the  x-axis  of  the  coordinate 
system  and  ZRUT  is  negative  to  transform  the  coordinates  of 
Z  to  their  zero  tilt  origin.  This  is  done  by 


2i  -  2i  ^ 


=  Ze 


-i$ 


or 


ZWORKB  =  ZWORKB*  ZRUT  . 

With  the  line  of  sight  defined  in  target  coordinates, 
the  azimuth  and  elevation  can  be  computed  as 

AZ  =  ATAN2 (WORKBR, TEMPO 
EL  =  ASIN (WORKBI*RRECIP) 

where  WORKBR  is  the  y  coordinate,  TEMPC  the  x  coordinate, 
WORKS I  the  z  coordinate  and  RRECIP  is  the  reciprocal  of  the 
range  magnitude. 
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Signal  Strength  Interpolation 

The  infrared  signal  strength  is  obtained  by  interpolat¬ 
ing  the  values  obtained  from  the  data  matrices  AZARAY  and 
ELARAY.  The  interpolation  is  done  by  a  comparison  of  values 
based  on  range#  azimuth  and  elevation  from  the  target.  Once 
a  representative  signal  is  calculated  it  is  compared  to 
minimum  standard  values  characteristic  of  the  missile  seeker 
head  to  determine  the  maximum  line-of-sight  tracking  rate. 

Azimuth 

The  program  initially  converts  azimuth  (AZ)  to  degrees# 
a  more  convenient  working  form  utilizing  the  radians  to 
degrees  conversion  factor  (RADDEG) .  An  index  (ITEMPA)  is 
next  computed  by  dividing  the  value  of  azimuth  by  15  and 
adding  one  to  it  so  the  lower  value  of  the  azimuth  band  can 
be  determined.  (The  azimuth  bands  run  from  the  nose  to  the 
tail  (0°  to  180°)  in  15°  bandwidths.) 

Next  an  increment  term  is  set  to  one  (IN  =  1)  for 
later  use.  If  the  azimuth  is  equal  to  180°  then  this  incre¬ 
ment  is  set  to  zero. 

TEMPB  (an  interpolation  angle  based  upon  the  size  of 
the  angle  formed  by  the  azimuth  with  the  lower  azimuth  band 
angle) ,  is  now  computed.  The  difference  between  the  azimuth 
and  the  lower  azimuth  band  angle  is  calculated#  and  this 
difference  is  multiplied  by  12  to  place  it  into  a  convenient 
size.  The  band  is  15°  wide  and  if  the  azimuth  equals  the 
lower  band  angle#  TEMPB  will  be  zero;  if  the  difference  is 
7.5°#  TEMPB  will  be  90°;  and  if  the  difference  is  15° # 
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TEMPB  will  be  180°.  TEMPA  is  then  set  to  the  cosine  of 
TEMPB  times  one  half. 

Data  points  from  the  azimuth  array  are  now  selected 
for  an  azimuth  interpolation.  TEMPC  is  chosen  by  the  com¬ 
puted  IRNG  and  ITEMPA  where  ITEMPA  is  the  lower  azimuth 
band.  TEMPD  next  selects  the  next  higher  azimuth  value  at 
the  same  range  as  TEMPC.  If  the  azimuth  is  180°,  IN,  which 
increments  ITEMPA,  is  zero  so  that  an  inappropriate  value 
is  not  selected  from  the  data  matrix.  (There  is  no  data 
point  for  azimuth  greater  than  180°.) 

The  interpolation  now  takes  place  at  the  lower  range 
band  and  a  lower  azimuth  band  value  (BNDAL)  is  computed: 

BNDAL  =  (TEMPC  +  TEMPD) * HALF  +  (TEMPC  -  TEMPD) * TEMPA 
or  equivalently: 

AZ_  Value  +  AZ.  Value 
BNDAL  =  — - - ^ - - - 

AZ  Value  -  AZ  Value 
+  — = - 2 - - -  cos  (TEMPB)  . 

If  TEMPB  is  0°  then  cosine (TEMPB)  is  equal  to  one  and  BNDAL 
reduces  to  the  lower  azimuth  value.  If  TEMPB  is  90°  then 
cosine  is  equal  to  zero  and  BNDAL  becomes  one-half  the  lower 
plus  one-half  the  upper  azimuth  values.  If  TEMPB  is  180° 
the  cosine  is  equal  to  (-1)  and  BNDAL  reduces  to  the  upper 
azimuth  value.  For  angles  between  these  values,  an  appro¬ 
priate  ratio  of  the  azimuth  values  is  computed. 

Next  an  upper  limit  (BNDAU)  based  on  the  next  range 
data  (IRNG  +1)  is  computed  in  a  similar  manner.  The  same 


azimuth  indices  are  used  but  IRNG  is  incremented  by  one 
into  the  next  outer  range  bin. 

Elevation 

The  next  series  of  computations  are  to  arrive  at  a 
pair  of  elevation  data  points  which  correspond  to  the  azi¬ 
muth  position.  In  other  words,  the  IR  signature  can  be 
thought  to  be  similar  along  a  conical  path  around  the 
source . 

The  first  statement  determines  if  the  elevation  is 
negative  and  if  it  ic#  adds  13  to  the  ITEMPA  index  to  ref¬ 
erence  the  negative  elevation  data  array.  Values  for  TEMPC 
and  TEMPD  are  selected  as  previously  discussed  for  azimuth, 
and  a  similar  interpolation  is  done  for  an  inner  range  value 
BNDEL  and  an  outer  range  value  BNDEU . 

The  final  interpolation  based  upon  angles  is  accomp¬ 
lished  by  first  computing  an  index  based  upon  the  elevation 
angle  (TEMPA) .  The  value  of  the  elevation  is  divided  by 
15  and  one  is  added  to  it  so  a  lower  value  of  the  elevation 
band  can  be  determined.  (The  elevation  bands  run  from  the 
plane  of  the  wings  to  the  perpendicular  axis  of  the  plane 
(0°  to  90°)  in  15°  bandwidths.) 

TEMPB,  an  interpolation  angle,  is  now  computed,  similar 
to  the  azimuth  interpolation  angle.  The  cosine  of  the  angle 
(TEMPB)  is  then  computed  and  the  final  interpolation  between 
the  azimuth  and  elevation  values  takes  place.  The  values 
computed  are  a  lower  (BNDL)  and  an  upper  (BNDU)  signal 
strength. 
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Final  Interpolation 

The  last  interpolation  combines  the  signal  values  at 
distances  which  bracket  the  actual  distance  at  the  correct 
azimuth  and  elevation.  The  missile  seeker  will  receive  a 
signal  somewhere  in  between  the  values  calculated  at  BNDL 
and  BNDU.  The  infrared  signal  is  a  function  of  one  over 
range  squared  (Ref  6) ,  so  what  the  program  does  now  is  com¬ 
pute  a  reference  signal  at  the  close  range  (BNDL)  and  then 
a  reference  signal  at  the  far  range  (BNDU) .  The  two  refer¬ 
ence  signal  strengths  are  then  combined  to  form  a  new  refer 
ence  signal  strength  at  the  actual  range,  and  finally  this 
reference  is  divided  by  the  range  squared  (TEMPD)  to  give 
the  interpolated  signal.  Let 


BNDU 

BNDL 

Then 


3y  *  Upper  tabulated  value  (Interpolated) 
S  «  Lower  tabulated  value  (Interpolated) 


and 


where  SrU  =  the  source  reference  signal  upper  and  SrL  *  the 
source  reference  signal  lower. 


Sru  -  BNDU*RBNDSQ(IRNG  +  1) 

SrL  =  BNDL*RBNDSC!(IRNG) 

Using  the  calculated  reference  strengths,  and  the  interpo¬ 
lated  signal  strengths  a  delta  signal  is  calculated  as  a 


79 


ratio  of  the  change  in  reference  to  the  change  in  interpo¬ 
lated  signals. 


AS 


(S 


rL 


£r0> 


<SU  '  SL> 


or 

TEMPB  =  (3NDL*TEMFC  -  BNDU*RBNDSQ ( IRNG  +  1) ) /BNDU  -  BNDL) 


TEMPA  is  now  computed: 

TEMPA  =  BNDL* (TEMPC  +  TEMPB) 
which  is  equivalent  to: 


TEMPA  =  SL*  (R^  +  AS)  . 

The  signal  is  then  calculated  by: 

RADANC  =  TEMPA/ (TEMPO  +  TEMPB) 


or. 


RADANC 


sJ^) 

L\rZ  +  AS/ 


which  is  the  value  of  the  signal  at  the  inner  range  decreased 
by  the  slope  determined  by  the  ratio  of  the  ranges  plus  the 
decreased  value  of  the  signal  from  the  inner  band  to  the 
outer  band. 

The  interpolated  signal  strength  is  then  converted  to 
Phasey  decibles  by: 


TENLNH  =  10 . G*ALOG (RADANC* PHASEY/TEMPD)  . 


where  Phasey  decibles  are  ten  times  the  natural  log  of  the 
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signal  and  PHASEY/TEMPD  is  a  factor  to  convert  the  interpo- 
lated  RAD AN C  from  watts/steradian  to  pico-watts/cm  .  The 
comparison  logic  uses  the  Phasey  decible  strength  to  deter¬ 
mine  the  maximum  track  rate  of  the  missile  seeker. 

Example  Problem 

For  an  example  of  how  the  azimuth  and  elevation  inter¬ 
polation  works,  suppose  the  target  and  launch  aircraft  have  a 
range  of  4000  feet  (IRNG=2) .  Further,  let  the  aspect  meas¬ 
ured  from  the  target's  nose  be  150°  and  the  elevation  be 
-45° . 

The  first  value  determined  is  ITEMPA  which  would  be 
equal  to  150  divided  by  15  or  10  plus  one  or  ITEMPA  would 
equal  11.  Since  150°  is  less  than  180°,  IN  would  equal  one. 

The  interpolation  angle  (TEMPB)  is  now  computed  by: 

TEMPB  =  (TEMP A  -  BNDAZ (ITEMPA) ) *12* DEGRAD  , 


or 


TEMPB  =  (150  -  150) *12*DEGRAD  . 

(DEGRAD  is  the  transformation  constant  from  degrees  to 
radians  just  as  RADDEG  is  the  transformation  constant  from 
radians  to  degrees.)  This  is  equivalent  to: 

TEMPB  =  0  . 

Next  TEMPA  is  computed  as  the  cosine  of  TEMPB  or 
TEMPA  =  COS (0) * { . 5)  =  0.5 
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New  the  first  interpolation  begins „  TEMPC  and  TEMPD  are 
the  values  which  bracket  the  azimuth  or  AZARAY  (XRNG, ITEMPA) 
and  AZARAY  (XRNG, ITEMPA  +  1).  Those  correspond  to  the 
values  165  and  180.  BNDAL  is  now  computed  by: 


BNDAL 


(TEMPC  +  TEMPD) 
2 


4* 


TEMPC  -  TEMPD 
2 


cos  (0°) 


or 


BNDAL  = 


165  +  180  +  165_ 


180  _ 


165 


BNDAU  is  computed  similarly  but  at  the  next  outer  range 
values  (IRNG  +  1) . 


BNDAU  = 


TEMPC  +  TEMPD  ,  TEMPC  -  TEMPD 


where  TEMPC  =  AZARAY  (3,11)  =  100  and  TEMPD  =  AZARAY  (3,12) 
=  110 .  So 


BNDAU  = 


100  +  110  +  100  ~  110 


=  100 


The  elevation  values  are  now  computed  in  a  similar 
manner,  but  the  elevation  is  less  than  0°  (-45°)  so. 


ITEMPA  =  ITEMPA  +  13 


or 


ITEMPA  =  24  . 

The  lower  elevation  value  is: 
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BNDEL 


ELARAY (2,24)  +  ELARAY (2,25) 

2  "" 


ELARAY (2,24)  -  ELARAY (2, 25) 


or 


BNDEL  = 


157  +  168  L  157  -  168 

^ 


=  157 


Similarly  for  the  outer  range: 


BNDEU  = 


105  +  109  L  105  -  109 

A  *  «  -  --  r-  ■  _ 


=  105 


The  final  computation  for  BNDL  and  BNDU  consists  of  first 
computing  an  interpolation  angle  as  a  function  of  the  ele¬ 
vation.  The  index  ITEMPA  is  computed: 

ITEMPA  =  (45/15)  +1  =  4. 

The  interpolation  angle  TEMPB  is  now  computed: 

TEMPB  =  (45  -  45) *12*DEGRAD 

and  the  factor  TEMPA: 

TEMPA  =  cos (0) * ( . 5)  =  0.5  . 

BNDL  is  now  computed: 


BNDL  = 


165  +  157  ,  165  -  157 

- - -  +  - s - 


=  165 


and 


BNDU  = 


100  +  105  ,  100  -  105 
- - -  +  - =— - 


=  100 
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Appendix  E 
DRIVER 


The  model  is  an  adaptation  of  a  dynamic  simulation 
utilized  by  the  ACMR/I  range  which  involved  real  time  data 
transmission  to  and  from  the  participating  aircraft  and  the 
simulation  computers.  To  enable  the  model  to  work  inter¬ 
actively  on  a  time  sharing  terminal,  a  DRIVER  program  had 
to  be  developed  which  would  initialize  the  necessary  param¬ 
eters  normally  supplied  by  the  aircraft  telemetry  pods.  In 
addition,  DRIVER  also  sets  up  the  quaternions  relating  the 
U  (earth  surface)  coordinate  system  to  the  T  (target) 
coordinate  system. 

Tilt  from  U  to  T 

The  tilt  quaternion  is  constructed  in  Section  2  through 
the  work  variables  TEMPA;  TEMPB;  WORKAR;  and  WORKAI,  as 
follows : 


TEMPA  =  VT  +  VTU1 
TUTR  =  VTUR/TEMPA 
TUTI  =  VTUI/TEMPA 

The  variables  used  here  are  VT,  the  magnitude  of  the  target 
velocity?  VTU1,  the  x-coordinate  of  target  velocity;  VTUR, 
the  z-coordinate  of  target  velocity;  and  VTUI,  the  y-coordi- 
nate  of  the  target  velocity.  These  lines  of  code  are  equiv¬ 
alent  to: 
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TEMPA  - 

VT 

+  VTUl 

VTUl  ^ 

vi 

t3 

-  TUTR 

= 

V2 

VT  +  V. 

VTUR  = 

V2 

fc2 

=  TUTI 

S3 

v3 

VT  +  v1 

VTUl  = 

v3  . 

Since  the 

tilt. 

t 

=  (1  =  kT) 

r 

and 


where 


so 


T  = 


t3  +  it2 


t„  = 


v. 


VT  +  v. 


and 


v. 


VT  +  v. 


T  =  TUTR  +  TUT I  =  ZTUT  . 


Roll  from  U  to  T 

The  roll  relation  is  formed  by  the  following  lines  of 


code: 


TEMPA  =  TUTI**  2 

DUT  =  TWO/ (ONE  +  TUTR* *2  +  TEMPA) 
WORKAR  =  ONE  -  DUT*TEMPA 
WORKAI  =  -DUT*TUTI*TUTR 
TEMPA  =  CABS  (ZWORKA) 

RUTR  =  -WORKAR/TEMPA 
RUTI  =  -WO  RKA I /TEMPA 

TEMPA  is  defined: 


2 


TEMPA 


(VT  +  v1) 


The  magnitude  of  T  is: 


(t2  +  t^) 1/2 


and  DUT  is  therefore: 


d  =  DUT  = 


1  +  |T| 


2  2 
1  +  4  +  *3 


1  +  (TUTI)2  +  (TUTR) 2 


The  additional  lines  of  code  follow  from  the  roll  transfor¬ 
mation 


M  -  A 


from  Appendix  B 


The  complex  number  A: 


A  =  1  -  i  d  Im(T)T 


where  again 

T 

f 

and  d  =  DUT 

ZRUT  which  is  equivalent  to  e1^  is  found  by: 

WORKAR  s  Re  (A)  =  1  -  dt2 

WORKAI  =  Im (A)  =  -dt2t3 
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and 


A  ~  ZWORKA  ■  WORKAR  4  iWORKAI  . 


Therefore 


A  -  1  -  dt^  -  idt2t3 


=  1  -  d (tf  +  it2t3) 


where 


i!m(T)T  = 


it2<t3  ' 


^2  +  ^2^3 


Since 


and 


TEMPA  =  CABS  (ZWORKA)  =  A 


RUTR  = 


RUTI  = 


Re  (A) 

WORKAR 

|A| 

TEMPA 

Im  (A) 

WGRKAI 

Ia! 

TEMPA 

then 


ZRUT  =  RUTR  +  RUTI  = 


A 

A 


Tilt  Derivative 

In  Section  three  of  DRIVER,  a  derivative  of  tilt  from 
U  coordinates  to  T  coordinates  is  calculated  based  upon  a 
roll  angle  <j>  from  the  target  to  the  maneuvering  plane?  NT, 
the  magnitude  of  the  normal  acceleration?  and  VT,  the 


SB 


magnitude  of  the  target  velocity.  The  code  begins  by  rede¬ 
fining  TEMPA: 

TEMPA  ®  TUTI**2 
which  is  equivalent  to: 

TEMPA  =  t\  . 


DUT  is  then  calculated: 

BUT  =  TWO/ (ONE  +  TUTR*  *  2  +  TEMPA) 


1  +  i  T  |  2 

WORKAR  ■  ONE  -  DUT* TEMPA 
WORKAI  =  -DUT*TUTI*TUTR 


d  =  DUT  =  - 5 - y 

1  +  t|  +  t2 

WORKAR  and  WORKAI  are  recalculated: 


or 

WORKAR  *  1  -  dt2 


and 


WORKAI  *  “dt2t3  * 
This  again  is  equivalent  to 

A  =  1  -  i  d  Im(T)T 
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or 


ZWORKA 

The  roll  relation  from 

ZRUT  = 


=  WORKAR  +  iWORKAI 

=  1  -  dtjj  “  idt2t2 

U  to  T  is  therefore: 

RUTR  +  RUTI  = 

|A| 

ZWORKA 
CABS (ZWORKA) 


The  derivative  of  the  tilt  is  now  calculated: 

TUTOR  =  (CPHITM*RUTR  -  SPHITM*RUTI) * TEMP A 
TUTDI  =  (-SPHITM*RUTR  -  CPHITM*RUTI) *TEMPA 

This  sequence  can  be  broken  up  into 

CPHITM*RUTR  -  SPHITM*RUTI 

where 

CPHITM  *  cost*) 

and 

SPHITM  »  sin4> 

<S>  being  the  angle  between  the  target  coordinate  plane  and  the 
maneuver  plane.  This  portion  of  the  code  projects  the  tar¬ 
get  roll  relation  onto  the  maneuvering  coordinate  frame  where 

CPHITM*RUTR  -  SPHITM*RUTI 
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is  the  real  axis  projection  and 


t 


-SPHITM*RUTR  -  CPHITM*RUTI 


is  the  imaginary  axis  projection.  The  sum  of  the  two  pro¬ 
jections  form  the  complex  number  O'  from  equation  (19)  in 
Appendix  B.  The  tilt  derivative  is  formed  by: 


*  = 


O' 

T 


This  is  accomplished  by  multiplying  the  terms  of  O'  by 
TEMPA  where: 


TEMPA  =  -NT/ (VT*DUT) 

NT 

where  the  factor  (-  — )  forms  the  derivative  of  O’  with 
respect  to  time  for  numerical  integration  and 

OUT  «  d 

which  completes  the  denominator  of  the  tilt  derivative: 

A  _  0^  _  TUTDR  -4-  TUTDI  ,  NT. 

d  DOT  VTJ 
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